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The exact computation of the correlation functions 
of ID quantum integrable models has been one of the 
challenging problems. For the spin-1/2 XXZ chain, the 
correlation functions in the antiferromagnetic regime 
were found to be expressed in the multiple integral 
form.^'^) However, the exact evaluation of them is still 
a hard work, and has been only successful for some spe- 
cial anisotropy parameters.'^-' 

In the ferromagnetic regime, there is a class of non- 
translationally invariant ground state which should be 
called as kink ground state**'^). We have studied the cor- 
relation function of the kink ground state, and exactly 
calculated the density matrix^^ (see also ref 7) . 

The kink ground state also exists for arbitrary spin 
and dimension. In this paper, we extend the analysis de- 
veloped in our previous paper to higher spin ID XXZ 
chain, and calculate the density matrix (note that the 
model we consider in this paper is not the integrable 
higher spin XXZ chain^) . 

The Hamiltonian of the spin S* ID infinite XXZ chain 

is 



H 



E 

mez 



qv qV 



(1) 

We shall consider the ferromagnetic regime A > 1. For 
later convenience, we parametrize A as A = {q^ + 
q~^)/2. Then A > 1 corresponds to < g < 1. A kink 
ground state is the superposition of kinks which have 
the same center. For a (normalized) kink state ^x&zlmx) 
(rrix e {-S, • ■ • , S}), the center j - 1/2 {j e Z) is the 
position where 



E 

x<j 



S), 



x>j 



holds. Let us denote the kink ground state whose center 
is at J — ^ by l^'j), and introduce the generating function 
of 



= (8) ( E 



' c{mx) 2 \mx)) 



® ( E N'^^+^^)™''+''c(m^)^I"^y)), 

i/eZ>o my = -S 

(2) 

where c{Tn^) = {28)1/(3 - m.^)l{S + m^)l. is the 
coefficient of of the expansion of |5'(z)), i.e, 



|vl'(z)}=^z^|vl/,-}. 



(3) 



Let us focus on one of the kink ground states |4'o), and 
calculate the density matrix 

where Ex'' ''\mx) = Sm^,e^Wx)- position of the 

site where the operator Ep-'^' acts on, and is assumed 
to be Xj Xk for f k. We only consider the case 
when X]J=i(^j ~ ^j) = is satisfied: otherwise, the 
density matrix is zero. We calculate (\E'(z)|\E'(z)) and 

nj=i -^^j*^' to obtain the exact expression 
of the density matrix since 



j = l j = -c>o j = \ 



(5) 



J/' 



(6) 



holds, and (l-ol^o), (*o| n"=i ^^^f l*o) can be ex- 
tracted from (*(z)|^'(z)} and (*(2)| 11^=1 -S^^'' l^(^)>> 
which are easier to calculate. Let us first calculate 
l^{z)\^{z)). Using the Jacobi triplet product identity 



('?;g)oc(-a;g2;g,)^(-x ^g2:g)^= ^ 



x^q 



(7) 



we have 
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where w = z^ , (a; q)oo := 11^10(1 ^ °-^^) ^^"^ 
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Using eq. (7) and the following identity^^ 
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{^{z)\ n^^i E^] ' \^{z)) can be calculated as 
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(10) 



j=-(x> jiei;,j2ei;>o 



where Cj = q^+''' , [e'jCj] = {ej + e'j + 2S)/2 and 



2S 



■jfe+n-l 



From eqs. (8) and (10), one has 



(*o|*o) 



(9; 9) 



2S • 



(11) 



X E ^hBn, (12) 

jiez,j2ez>o 
ji+j2=-E7=iK-'^j] 

leading to the exact expression of the density matrix 
3=1 j=i 

E (13) 

jiez,j2ez>o 

jl+j2 = -E"=l[4«j] 

One can check that for S = 1/2, cq. (13) recovers the 
result in ref. 6. However, the expression of eq. (13) gets 
more complicated as the spin becomes higher. 

Concentrating on the S = 1 case, we can derive some 
slightly easier expression by using 



n ^ 00 

j=l ^ 3=0 



which is a special case of eq. (9). The result is 

(n4^^)=n#'^^^(^^)^(4)} 



3 = 1 



i=i 



00 



X ^(-l)-'X„j93a+EJ=i[4'=/=])' 

3=0 

x(5^"<=" .+06"%.^ r, 1), (14) 

where C = 2 ^'V(l + 2 E^i ?'"')• From cq. 

(14), the magnetization and the spin-spin correlation 
functions can be easily calculated. 

(^.^) = E(-i)^(j + i)C^a^ 
3=0 

X (C^g^+1 - l)((5f + C5f'^), (15) 

00 

(S'A) = ^(-l)^X2,,(5f- + CSf^) 
i=o 

X {CiCiq— - id + CI)q~), (16) 



(5+S-)=4Ci^C2^E(-l)'^2,,- 



j=0 



(3 + 2)^ 



X {(^f'=" + C5f'')(Ci +C2)g 

+ (,5^^ + CJ|-")(g^ + CiC2g^)}, 



(17) 



where 



_ (i + i)(cr' - CVl - (j + 3)CiC2(Cf^^ - CD 
(Ci-C2)3 

From these expressions, one can easily show that the 
spin-spin correlation functions decay exponentially for 
large distances. 

{S'S^.^ - (5|,)(5^2) ~ A^'{xi)q^'+i for X2 » 1, 



(18) 



^-(a;i)=2^(-l)^9^(l-CiV+') 

3=0 

X (jCr' + U + i)g^CCf)(<5f ^" + C(5f '^), 

(5+ S- ) ~ A+-(xi)<7^("^+^) for » 1, 

00 

^+-(a;i) = 4Cf ^(-Ci)^'(i + 1) 

j=0 



(19) 
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